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ABSTRACT 


A  nonlinear  maneuvering  model  based  on  ship  geometric  and  mass  properties  is 
developed.  The  model  can  be  utilized  to  evaluate  maneuvering  performance  early  in 
the  design  phase.  This  model  is  also  used  in  this  thesis  as  a  benchmark  in  order  to 
evaluate  the  accuracy  of  the  simpler  Nomoto's  model.  The  latter  is  faster  to  simulate 
and  is  ideally  suited  for  visual  simulation  studies.  Results  comparing  the  relative 
accuracy  and  speed  of  implementation  of  the  two  models  are  presented  for  different 
inputs  and  geometric  properties.  An  improvement  to  Nomoto's  model  is  suggested 
which  greatly  increases  accuracy  while  maintaining  high  speed  of  real  time 
implementation.  Furthermore,  a  series  of  parametric  studies  is  performed  in  order  to 
evaluate  the  sensitivity  of  fundamental  maneuvering  properties  in  terms  of  basic  ship 
geometric  quantities. 
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I.  INTRODUCTION 


Today's  military  faces  many  challenges  to  do  more  with  less.  Since  the  end  of  the 
Cold  War,  the  topic  of  discussion  on  Capitol  Hill  has  been  the  shrinking  military  budget. 
However,  the  military's  requirement  to  maintain  full  combat  readiness  has  not  changed.  In 
order  to  maintain  training  at  the  required  level,  alternatives  to  costly  full  scale  maneuvers 
must  be  developed. 

Current  war  gaming  models  use  mostly  static  or  simplified  dynamic  descriptions  of 
both  attacking  weapons  and  targets.  Clearly,  this  provides  an  unrealistic  description  of  the 
situation,  with  respect  to  the  vulnerability  of  the  target  to  attack,  and  does  not  allow  for  the 
incorporation  of  standard  tactical  maneuvers  of  the  target  to  avoid  incoming  weapons.  What 
needs  to  be  done  is  to  increase  the  accuracy  of  current  models  by  incorporation  of  targets  in 
a  dynamic,  rather  than  a  static  manner.  Developing  an  accurate  simulation  model,  will  lead 
to  ship  designers  being  able  to  examine  the  vulnerability  of  a  ship,  in  real  time  through  the 
simulation  model. 

In  this  thesis,  a  nonlinear  model  based  on  ship  geometric  and  mass  properties  will  be 
developed.  The  model  can  be  utilized  to  evaluate  maneuvering  performance  early  in  the 
design  phase.  This  model  will  also  serve  as  a  benchmark  in  order  to  evaluate  the  relative 
accuracy  of  a  simpler  linear  model.  Chapter  n  will  show  the  development  of  the 
mathematical  models.  Chapter  HI  will  discuss  the  assessment  of  the  models  and  compare 
them  for  visual  simulation  suitability.  Chapter  IV  will  evaluate  the  sensitivity  of 
fundamental  maneuvering  properties  in  terms  of  basic  ship  geometric  quantities.  Chapter 
V  will  assess  two  variations  to  the  linear  model,  and  compare  them  to  the  nonlinear  model. 
Finally,  Chapter  VI  will  discuss  conclusions  and  recommendations  for  further  work. 
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II.  MATHEMATICAL  MODELING 


In  this  dynamical  model,  the  planar  motion  of  a  ship  on  an  undisturbed  free-suiface 
of  an  incompressible  fluid  will  be  considered.  The  basic  dynamics  of  maneuvering  will  be 
described  and  analyzed  using  Newton's  equations  of  motion. 

We  will  define  two  reference  frames,  one  set  of  axes  fixed  relative  to  the  Earth 
aligned  with  directions  North,  East,  and  Down  (denoted  by  capital  letters),  and  a  second  set 
of  axes  fixed  relative  to  the  ship,  aligned  with  the  longitudinal,  the  starboard  direction,  and 
down  (denoted  by  lower  case  letters). 


Figure  1  Designation  of  reference  frames 


The  heading  (yaw)  angle,  refers  to  the  direction  of  the  ship's  longitudinal  axis,  with 
respect  to  one  of  the  fixed  axes.  The  drift  angle  refers  to  the  difference  between  the  heading 
and  the  actual  course.  Let  the  instantaneous  speed  of  the  ship  be  denoted  by  U,  then  the 
surge  and  sway  velocities  are  given  by: 

u^C/cosp  ,  v=-t/si«p 
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The  dynamical  model  used  in  this  work  has  been  proposed  [Ref.  1]  and  is  given  by; 


Y'—(m  '^m^L/UXUsin^/U*Pco&^y{m  '*fn^  'cosp 


Ar'=  ii;,.i;,)(L/U\Ur\Ur) 


where, 

P  is  the  drift  angle 
r  is  the  rotational  or  yaw  velocity 

,and  i^  are  the  mass  moments  of  inertia  in  both  ref  frames 
m,  m,(,  and  niy  are  the  ship's  mass,  added  mass  of  the  x-axis,  and  the 
added  mass  of  the  y-axis,  respectively 

AH  terms  of  these  equations,  and  the  equations  of  motion  in  the  remaining  report,  have  been 
nondimensionalized  using  the  ship's  length  L,  the  ship's  draft  d,  the  fluids  density  p,  and  the 
sh^'s  speed  U,  as  diown  in  the  Appendix  The  terms  X,  Y,  and  N  are  the  external  force  along 
the  x-axis,  the  external  force  along  the  y-axis,  and  the  yaw  moment  about  the  ship's  center  of 
gravity,  respectively. 

It  is  assumed  that  X,  Y,  and  N  have  the  following  decomposition, 

X.X^Xp.X^ 

^•Yif*Tp*Tp 

N-Njf*Np^Np 

where  the  subscripts  represent  respective  contribution  due  to, 

H  -  ship's  hull 
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P  -  ship's  propeller 
R  -  ship's  rudder 

We  can  assume  that  during  most  ordinary  maneuvers,  the  forward  velocity  is  kept  constant 
by  the  propulsion  control.  Also,  in  looking  at  the  ship's  response  to  rudder  deflection,  it  is 
assumed  that  no  significant  surge  accelerations  will  take  place.  From  these  assumptions  we 
are  left  with, 

N-N„^N, 


The  components  of  the  rudder  forces  are  taken  to  be  of  the  form. 


F^=-(^«/^Wsin(6-(P-2x^)) 


C^=6.13^/(is:^.2.25) 


where. 


F^-  force  acting  on  rudder 

Ar-  rudder  area 

Kr-  rudder  aspect  ratio 

Ur-  effective  rudder  inflow  speed  (taken  to  be  equal  to  U) 

Xr-  distance  between  COG  of  ship  and  center  of  lateral  force 
Xh-  distance  between  COG  of  ship  and  center  of  additional  lateral 
force 

aH*  ratio  of  additional  lateral  force 
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Ar.  Kr,  Xr,  Xh,  and  an  are  constants  obtained  through  model  tests. 


The  lateral  force  and  yaw  moment  acting  on  the  ship’s  hull  are  expressed  as  follows: 


The  lateral  force  and  yaw  moment  hydrodynamic  coefficients  are  given  as, 

r=m^-l.5sCg 

r.=-TzkAAsC„ 

p  2  B 

V2.97i(l-Q 


r,p=1.25^(l-C>0.5 

Y  =0.17-C„-0.07 
7p,^=0.75^C„-0.65 


N^-k 

JVpp=-0.48A(l-C>0.66 


N  =0.55C„-0.09 

rr  a 


A^^„=-0.25^C,.0.05 

Arpp,=-57.5(^C,M8.45C3-1.6 
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where, 


s=BIL,  Slenderness  ratio 


k=2dlL,  Aspect  ratio 


Block  coefficient 


m^=^2sCJi\*k),  Mass  coefficient 


The  above  model  may  be  simplified  by  linearizing  the  sway  and  yaw  equations,  as  done  in 
[Ref.  3],  to  obtain. 


iI^-Ny-(N,-mx^v-.N^v^iN^-mx^U)r^N,b 


where  6  is  the  rudder  deflection  angle  measured  according  to  the  convention  of  Figure  (1); 
positive  rudder  corresponds  to  a  turn  to  port.  These  equations  of  motion  can  be  manipulated 
to  produce  a  pair  of  decoupled  second  order  equations  in  r  and  v. 


TJ^HTcT,)v*v-Kb*KTfi 


We  will  assume  that  there  is  no  appreciable  side  slip,  and  hence,  the  second  equation  may 
be  ignored.  What  remains  is  known  as  Nomoto's  equation.  The  two  indicies,  K  and  T,  are 
related  to  the  hydrodynamic  coefficients  by  the  following, 

^  iY^-mXK-mx^U)<N,-IJY^ 

'■  YSN-mx^U)^pr-mU) 
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where, 


Y^(N-mx^U)-NXY^-mU) 

^  (N^-mx^Y,-(Y^-m)N, 

jjr_ _ V  a  V  6 _ 

~  YXN-mx^U)-N(Y^-mU) 

Y^^N^O.O 

r=m,-1.5.C, 

r=-0.5n/t4.45C3 

Y,Hya„)F^ 

K=-k 
NJ  -I 

f  u 

m^^2sC^(Uk),  i-\,2  k^=t,  A~(1-4); 

s  2 
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o  s 

Nomoto's  equation  expresses  the  relationship  between  the  ship's  turning  rate  and  the  rudder 
angle  through  a  second  order  transfer  function. 


K^KT^s 

which  can  be  simplified  to  a  first  order  transfer  function, 

r  K 
6'TsA 

This  equation  in  differential  form  is  known  as  Nomoto's  First  Order  Equation, 

mr-Kb 


where. 


By  using  basic  principles  of  ordinary  differential  equations,  the  solution  for  r,  in  terms  of  K 
and  T,  is  given  by. 


r  =  k:6(l-e  "O 


Having  developed  these  two  models,  nonlinear  and  liner,  the  two  will  now  be  compared  to 
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see  which  will  give  a  faster  calculation  time,  and  which  is  more  accurate  for  graphical 
simulation. 
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III.  MODEL  ASSESSMENT 


Having  developed  the  nonlinear  and  linear  (Nomoto's)  models,  the  two  were 
compared  for  speed  of  implementation  and  the  relative  accuracy  for  different  inputs  and 
geometric  properties. 

The  principal  dimensions  of  length,  beam,  and  draft  were  varied  independently  and 
the  models  were  also  compared  utilizing  different  rudder  angle  inputs. 

In  every  instance  of  the  comparison  between  these  two  models,  the  linear  model 
proved  to  be  approximately  six  times  faster  than  the  nonlinear  model.  This  fact  makes 
Nomoto's  model  ideally  suited  for  visual  simulation  studies.  Visual  simulations  may  be 
carried  out  with  minimal  threat  of  reduced  frame  speed,  due  to  calculation  processing  time. 

A  graph  of  the  advance  vs.  transfer  (i.e.,  a  geographical  x-y  plot  of  the  ship's  path) 
is  shown  in  Figure  2  for  both  linear  and  nonlinear  models  and  for  different  rudder  angles. 
Advance  and  transfer  are  shown  in  dimensionless  form,  with  respect  to  the  ship  length. 
Rudder  angle,  indicated  by  d  on  the  graphs,  is  in  degrees.  Solid  curves  correspond  to  the 
nonlinear  simulation,  while  dotted  curves  correspond  to  Nomoto's  first  order  model.  It  can 
be  seen  that  the  latter  accurately  depicts  the  ship's  path  for  relatively  large  rudder  angles. 
However,  as  the  rudder  angle  is  decreased  the  difference  between  linear  and  nonlinear 
simulations  is  likewise  more  pronounced.  This  deviation  may  be  attributed  to  the  increased 
side-slip  velocities  that  are  encountered  during  smaller  rudder  deflections  and  are  neglected 
in  Nomoto's  model.  This  is  explored  in  the  next  chapter. 

The  results  for  the  yaw  rate  r  are  shown  in  Figure  3.  Both  yaw  rate  and  time  are 
shown  in  dimensionless  quantities,  with  respect  to  ship  speed  and  length.  It  seems  that 
higher  rudder  angles  result  in  slight  overshoot  in  the  yaw  rate  for  the  nonlinear  model.  This 
means  that  for  high  rudder  angles,  the  yaw  rate  is  basically  described  by  a  second  order 
system,  whose  damping  ratio  becomes  smzdler. 
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The  sensitivity  of  these  results  for  small  variations  in  the  ship's  draft  is  demonstrated  by 
Figures  4  through  7.  A  plot  of  the  advance  vs.  transfer  for  different  values  of  the  draft  based 
on  the  nonlinear  model  is  shown  in  Figure  4.  It  appears  that  small  drafts  correspond  to  more 
maneuverable  ships  with  tighter  turning  circles.  The  same  result  can  be  drawn  from  the 
linear  simulation  results  of  Figure  5.  From  the  yaw  rate  simulations  of  Figures  6  and  7,  we 
can  see  that  smaller  drafts  correspond  to  higher  turning  rates  at  steady  state,  as  well  as,  at 
smaller  time  constants.  The  nonlinear  simulations  also  suggest  that  a  smaller  draft  results 
in  smaller  damping  ratios,  as  evidenced  by  the  overshoot  in  the  yaw  rate  time  history. 

The  sensitivity,  with  regards  to  ship  beam,  is  demonstrated  by  Figures  8  throug:  11. 
In  general,  it  appears  that  the  beam  has  very  little  effect  on  maneuverability.  The  nonlinear 
model  suggests  that  maneuverability  is  increasing  for  increasing  beam,  although  the  linear 
model  predicts  the  opposite.  It  also  appears  that  the  yaw  rate  damping  ratio  is  decreased  for 
decreasing  values  of  the  ship's  beam. 

The  sensitivity,  with  regards  to  the  ship  length,  is  demonstrated  by  Figures  12  through 
18.  In  general,  both  linear  and  nonlinear  models  suggest  that  maneuverability  is  increased 
with  increasing  length.  This  results  in  tighter  turning  circles,  both  in  dimensional  and 
dimensionless  quantities. 
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Figure  2  Rudder  Angle  comparison  between  Nonlinear  Model  and 
Nomoto's  Model 


Figure  3  Yaw  Rate  comparison  between  Nonlinear  Model  and  Nomoto's 
Model 
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Figure  ^  Advance  vs.  Transfer  for  Nonlinear  Model  at  different  drafts 


Figure  5  Advance  vs.  Transfer  for  Nomoto's  Model  at  different  drafts 
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Figure  8  Advance  vs.  Transfer  for  Nonlinear  Model  at  varying  beams 


Figure  9  Advance  vs.  Transfer  for  Nomoto's  Model  at  varying  beams 
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Figure  12  Advance  vs.  Transfer  for  Nonlinear  Model  at  varying  lengths 


Figure  13  Dimensional  Advance  vs.  Transfer  for  Nonlinear  Model  at 
varying  lengths 
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Figure  14  Advance  vs.  Transfer  for  Nomoto's  Model  at  varying  lengths 
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IV.  LINEAR  MODEL  VARIATION  ASSESSMENT 


In  Chapter  II,  the  assumption  of  negligible  side  slip  velocity  was  made  to  obtain 
Nomoto's  equation.  If  this  assumption  was  not  made,  we  would  have  the  pair  of  decoupled 
second  order  equations  from  the  linearized  sway  and  yaw  equations. 

The  relationship  between  the  indicies,  in  the  yaw  rate  second  order  equation  and  the 
hydrodynamic  coefficients  was  shown  in  Chapter  D.  The  indices  of  the  side  slip  velocity 
second  order  equation  are  related  to  the  hydrodynamic  coefficients  by  the  following, 

riN-mx^U)-NJir-mU) 

(N-mx^U)Y,-{Y-mU)N, 

The  second  decoupled  second  order  equation  led  to  formulation  of  Nomoto's  equation,  and 
the  first  expresses  the  relationship  between  the  side  slip  velocity  and  the  rudder  angle  throu^ 
a  second  order  transfer  function. 


which  can  be  simplified  to  a  first  order  transfer  function, 

V 

6’rs4 

V 
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This  equation  in  differential  form  reduces  to, 


Tv^v-^K6 

V  V 


where, 


Once  again,  by  using  basic  principles  of  ordinary  differential  equations  the  solution  for  the 
side  slip  velocity  is. 


v-A:^6(l-e  ) 

In  Figure  18,  a  comparison  is  shown  between  Nomoto's  model  and  a  linear  model 
containing  the  affects  of  side  slip.  The  two  graphs  are  almost  identical,  which  leads  to  the 
conclusion  that  the  original  assumption  of  negligible  side  slip  velocity  was  correct. 

Now  we  will  compare  the  nonlinear  model  with  the  side  slip  model.  Figures  19  and 
20  lead  to  the  same  conclusions  as  Figures  2  and  3,  that  relatively  speaking,  the  nonlinear 
model  is  more  accurate  than  the  linear  model.  Figure  21  shows  the  comparison  between  the 
side  slip  velocity  obtained  in  the  nonlinear  vs.  the  linear  case.  Once  again,  it  is  clear  that  the 
nonlinear  model  simulates  the  situation  more  accurately. 

The  side  slip  assessment  leads  to  the  conclusion  that  the  difference  between  the  linear 
and  nonlinear  models  lies  not  with  the  side  slip  velocity,  but  with  the  yaw  rate.  It  can  be  seen 
in  Figure  22  that  the  difference  in  the  error  between  the  linear  and  nonlinear  models  is  a 
function  of  the  rudder  angle. 

By  plotting  the  rudder  angle  vs.  the  steady-state  yaw  rate  ratio  (nonlinear  over  linear), 
we  obtain  a  straight  line  as  shown  in  Figure  23.  By  calculating  the  slope  and  intercept  of  this 
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line  and  applying  it  as  a  correction  factor,  we  obtain  a  new  relationship  for  Nomoto's 
equation, 

Tf  .  r=  (a-p|61)A:6 

By  applying  this  "correction"  factor,  it  can  be  seen  in  Figure  24  that  we  obtain  a  model  that 
is  more  accurate  than  both  Nomoto’s  model  and  the  side  slip  model.  On  the  average  the 
corrected  model  is  approximately  10%  slower  than  Nomoto's  model. 

In  Figure  25,  it  can  be  seen  that  the  relative  error  in  the  yaw  rate  between  the 
nonlinear  model  and  Nomoto's  model  has  been  removed  in  the  corrected  model.  By 
comparison,  the  conclusion  that  the  corrected  model  is  the  most  accurate  of  the  models 
presented  can  be  made.  Due  to  the  relative  accuracy  of  the  corrected  model,  it  is  better 
suited  to  simulate  the  turning  rate  while  maintaining  an  acceptable  frame  speed  during  visual 
simulations. 
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Figure  18  Advance  vs.  Transfer  comparison  between  Nomoto's  Model  and 
Side  Slip  Model 


Figure  19  Advance  vs. Transfer  comparison  between  Nonlinear  Model  and 
Side  Slip  Model 
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Figure  20  Yaw  Rate  vs.  Time  comparison  between  Nonlinear  Model  and 
Side  Slip  Model 


Figure  21  Side  Slip  Velocity  comparison  between  Nonlinear  Model  and 
Side  Slip  Model 
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-  NONLINEAR 


—  LINEAR 


Figure  22  Yaw  Rate  vs.  Time  comparison  between  Nonlinear  Model  and 
Nomoto's  Model  at  varying  rudder  angles 


Figure  23  Yaw  Rate  Ratio  vs.  Rudder  Angle  plot  used  to  obtain  correction 


factor  for  Nomoto's  Model 
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Figure  25  Yaw  Rate  vs.  Time  plot  comparing  Nonlinear,  Nomoto,  and 
corrected  Nomoto  Models 
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V.  PARAMETRIC  STUDY 


During  the  preliminary  ship  design  phase,  it  is  often  possible  to  introduce  small 
changes  in  fundamental  hull  properties  so  that  maneuverability  is  improved.  Of  course,  such 
changes  may  only  be  instituted  after  taking  due  consideration  of  all  aspects  and  requirements 
of  the  design  problem. 

In  this  Chapter,  the  maneuvering  properties  of  a  ship  are  explored  by  utilizing  the 
nonlinear  model  described  in  the  previous  chapters,  since  it  reflects  reality  better  than  the 
linear  model.  In  order  to  establish  a  systematic  series  of  runs,  the  ship's  displacement  was 
maintained  constant.  The  dispalcement  is  usually  determined  by  payload  considerations  and 
it  is  unlikely  that  a  designer  would  have  the  freedom  to  change  it  based  on  maneuvering 
requirements  alone. 

Figures  26  and  27  show  results  that  are  obtained  by  keeping  both  the  displacement 
and  the  beam  constant.  The  length  and  draft  are  then  changed  within  plus/minus  5%  of  their 
nominal  values.  It  can  be  seen  that  the  longer  and  shallower  hull  form  has  superior 
maneuvering  characteristics  compared  to  shorter  and  deeper  draft  ships. 

Figures  28  and  29  show  results  that  are  obtained  by  keeping  both  the  displacement 
and  the  length  constant.  The  beam  and  draft  are  then  changed  within  plus/minus  5%  of  their 
nominal  values.  It  can  be  seen  that  wide,  shallow  hull  forms  have  superior  maneuvering 
characteristics  compared  to  narrow  and  deeper  draft  ships. 

Figures  30  and  3 1  show  results  that  are  obtained  by  keeping  both  the  displacement 
and  the  draft  constant.  The  length  and  beam  are  then  changed  within  plus/minus  5%  of  their 
nominal  values.  It  can  be  seen  that  the  long  and  narrow  hull  forms  have  superior 
maneuvering  characteristics  compared  to  shorter  and  wider  draft  ships. 

Maneuvering  performance  is  increased  by: 

-  a  longer  length 

-  less  draft 

-  a  narrower  beam  for  constant  draft 

-  a  wider  beam  for  constant  length 
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Therefore,  beam  does  not  have  a  consistent  trend  and  since  it  is  usually  set  by  static 
stability  considerations,  it  should  not  be  altered  based  on  maneuvering  alone.  The  length  is 
the  most  costly  dimension  to  increase.  Not  only  are  material  cost  increased,  but  dry  dock 
size,  and  shipbuilding  considerations  must  also  come  into  play.  The  draft  on  the  othe  hand 
will  have  it's  own  set  of  considerations  such  as  harbor  and  canal  restrictions. 
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Figure  26  Advance  vs.  Transfer  comparison  with  constant  beam 


Figure  27  Yaw  Rate  vs.  Time  comparison  with  constant  beam 
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Figure  28  Advance  vs.  Transfer  comparison  with  constant  length 


Figure  29  Yaw  Rate  vs.  Time  comparison  with  constant  length 
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Figure  30  Advance  vs.  Transfer  comparison  with  constant  draft 


Figure  31  Yaw  Rate  vs.  Time  comparison  with  constant  draft 
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VI.  CONCLUSIONS 


In  this  thesis  two  maneuvering  models  were  developed  and  compared.  A  nonlinear 
model  based  on  ship  geometric  and  mass  properties  was  developed  from  [Ref.  1].  The 
second,  linear,  model  was  developed  by  linearizing  the  sway  and  yaw  equations  of  the 
nonlinear  model,  as  done  in  [Ref  3].  The  models  were  compared  for  speed  of 
implementation  and  the  relative  accuracy  for  different  inputs  and  geometric  properties.  It 
was  concluded  that  though  the  linear  model  was  approximately  six  times  faster  than  the 
nonlinear  model,  it  was  not  accurate  enough. 

At  this  point  a  variation  was  made  to  the  linear  model  by  adding  the  effects  of  the 
side  slip  velocity  to  possibly  explain  the  difference  in  the  accuracy.  Upon  further  inspection, 
it  was  determined  that  the  difference  in  accuracy  did  not  lie  with  the  side  slip  velocity,  but 
with  the  difference  in  the  yaw  rates. 

By  plotting  the  rudder  angle  vs.  the  steady-state  yaw  rate  ratio,  a  correction  factor  was 
obtained  and  applied  to  the  linear  model  ;causing  the  steady-state  error  between  the  yaw  rates 
of  the  linear  and  nonlinear  models  to  be  removed  altogether.  This  "corrected"  linear  model 
also  proved  to  be  relatively  more  accurate  than  either  the  original  linear  model  and  the 
nonlinear  model.  Though  this  model  is  10%  slower  than  the  original  linear  model,  it  is  still 
much  faster  than  the  nonlinear  model. 

A  series  of  parametric  studies  was  also  performed  on  the  nonlinear  model  in  order 
to  evaluate  the  sensitivity  of  fundamental  maneuvering  properties  in  terms  of  basic  ship 
geometric  quantities.  It  was  concluded  that  maneuvering  performance  is  increased  by  a 
longer  length  or  by  less  draft.  The  beam  did  not  show  a  consistent  trend,  and  is  usually  set 
by  static  stability  considerations  and  probably  not  alterable  for  maneuvering  alone. 

For  further  work,  it  is  recommended  that  these  models  be  incorporated  into  the 
Computer  Science  Department  battlespace  simulation  system,  NPSNET.  Initially  this  system 
was  geared  towards  ground  forces  and  land  based  conflicts,  but  recently  a  naval  component 
has  been  added  in  an  attempt  of  attaining  a  more  realistic  joint  services  simulator.  The 
models  should  be  checked  for  accuracy,  realism,  and  speed  of  implementation.  Once  a 
baseline  is  obtained,  improvements  in  the  best  model  for  visual  simulation  can  be  made. 
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APPENDIX.  MATLAB  CODE 


Nondimensional  quantities,  denoted  by  primes,  are  defined  as  follows: 


m'=ni/(0.5pLM) 

mx'=“x/(0-5pLM) 

my'=my/(0.5pL^d) 

IJ=I^/(0.5pLM) 

W=y(0.5pLM) 

X'=X/(0.5pLdU^) 

Y’=Y/(0.5pLdU') 

N’=N/(0.5pLMU^) 

r'=rL/U 
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%  COMPARISON  OF  LINEAR  AND  NONLINEAR  MODELS 

%  ON  SAME  GRAPH { including  Nomoto  correction  U  side  slip  model) 
% 

%  Define  parameters  of  the  simulation  (Metric  units) 

% 


time  =  90.0; 
dt  =  .1; 

Iz  =  3257400. ; 

L  =  250.0; 

RHO  =  1025; 

XG  =  10.0; 

Mass  =  147241915.0; 
m  =  Mass/ (0 . 5*RHO*L''3)  ; 

Izz  =  Iz/  (0 . 5*RHO*L'"5)  ; 

Xg  =  XG/L; 
d  =  -0.2618; 

B  =  40.77; 

D  =  20.00; 

S  =  B/L; 

V  =  Mass /RHO; 

Cb  =  V/ (L*D*B) ; 
k  =  2-^D/L; 

Kr  ::  1.7: 

AH  ^  74.0; 

D;<  ^  0.5: 

P  :  0  B*Dp; 
n  : :  j  0  0  0  ; 

Ctp  ^  1.  Oc¬ 
elli  =-•  0.63  3  *Cb  -  0.153; 

c;p.s  ^  -156  -  2*  (Cb*B/L)  ^2  +  4 1 . 6 *  ( Cb*B/L ) 
wp  --  0.23  +  1.4*  (CP  -  0.5)  "2; 
cp  -  0.6*wp; 
wr  =  1  -  (eps*  ( l-\Arp)  )  ; 

K  =  0.6*(l-wp)/(l-wr) ; 
tr  =  0.28*Cb  +  0.55; 
eta  =  Dp/ (sqrt (AR*Kr)  )  ; 

Cn  =  (6.13*Kr)/{Kr  +  2.25) ; 
lambda6  =  (k/s ) *  ( 1-1 . 6 *s ) ; 

I  =  (1/8) *s*Cb* {l+lambda6) ; 
lambda2  =  (k/s ) * ( l-s/2 )  ; 
m2  =  2*s*Cb* (l  +  lambda2)  ; 
lambdal  =  k; 

ml  =  2*s*Cb* ( 1  +  lambdal )  ; 

PI  =  3.14; 
xR  =  -0.5; 
xH  =  -0.5; 

ml  =  2 *s*Cb* ( 1  +  lambdal )  ; 

gamma  =  -22 . 2  * (Cb*B/L) ^2  +  0 , 02  * (Cb*B/L) 
U  =  1.0; 

UR  =  U; 


%  Final  time  of  simulation 
%  Sample  time  of  difference  eqns 
%  Moment  of  inertia 
%  ship's  length 
%  Density  of  salt  water 
%  Center  of  gravity 
%  Ship  mass 

%  Dimensionless  ship  mass 
%  Dimensionless  moment  of  inertia 
%  Dimensionless  center  of  gravity 
%  rudder  angle 

%  ship's  breath (twice  the  beam) 

%  ship's  draft 
%  slenderness  ratio 
%  ship's  displacement 
%  block  coeff 
%  hull  aspect  ratio 
%  rudders  aspect  ratio 
%  rudder  area 
%  propeller  diameter 

%  propeller  revolution  (rps) 


-  1.76; 

%  thrust  deduction  coefficient 


%  LINEAR  MODELS 
% 

beta  =  .25;  %  drift  angle 

Fn=  - (AR/ (L*D) ) *Cn*UR^2 *cos ( -  gamma*beta) ;  %  Force  acting  on  rudder 
% 

%  Define  hydrodynamic  coefficients 
% 

Yvl  =  -0.5*PI*k  +  1.4*s*Cb; 

Nvl  =  -k; 

Yrl  =  ml  -  1.5*s*Cb; 

Ydl  =  -(l+aH)*Fn; 

Nrl  =  -0.54*k  +  k"^2; 

Ndl  =  -(xR  +  aH*xH)*Fn; 

Nrl dot  =  Izz-I; 

Yvl dot  =  m-m2 ; 

Yrldot  =  0.0; 

Nvldot  =  0.0; 

% 

%  Define  constants  of  Nomoto 's  first  order  equation 
% 

K  =  (Nvl*Ydl  -  Yvl*Ndl) / [Yvl* (Nrl-m*Xg)  -  Nvl* (Yrl-m) ] ; 

Kv  =  - [ (Nrl-m*Xg) *Ydl  -  ( Yrl-m) *Ndl ]/ [Yvl * (Nrl -m*Xg )  -  Nvl* (Yrl-m) ] ; 

T1  =  [ (Yvldot-m) * (Nrl-m*Xg)  +  Yvl* (Nrldot-Izz)  ] / [Yvl* (Nrl-m*Xg)  -  Nvl* (Yrl-m)  ] 
T2  =  [  (Nvldot-m*Xg)  *  (Yrl-m)  -  Nvl*  (Yrldot-m*Xg)  ]  /  [Yvl*  (Nrl-m*Xg)  -  Nvl*  (Yrl-m) 
T3  =  [Ydl* (Nvldot-m*Xg)  -  Ndl* (Yvldot-m) ]/ [Nvl*Ydl  -  Yvl*Ndl]; 
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T^="TliT2-T3^  '  ■  (Yrldot  -  m*Xg) *Ndl ) / [ (Nrl-m*Xg) *Ydl  -  ( Yrl-m) ‘Ndl J ; 

Tv  =  T1+T2-T4; 

% 

%  Define  storage  vectors  for  Nomoto's  Model 
% 


kmaxl  =  (time/dt)+  1; 
rl  =  zeros ( 1 , kmaxl ) ; 
rldot  =  zeros  {1,  Janaxl)  ; 
psil  =  zeros  ( 1 , kmaxl ) ; 
psildot  =  zeros { 1 , kmaxl ) ; 
XL  =  zeros ( 1 , kmaxl ) ; 

XLdot  =  zeros  ( 1 , kmaxl ) ; 

P) 

YL  =  zeros ( 1 , kmaxl )  ; 

YLdot  =  zeros ( 1 , kmaxl ) ; 
lip) 

% 


%  Number  of  data  points 
%  Initialize  angular  velocity  vector 
%  Initialize  angular  acceleration  vector 
%  Initialize  heading  angle  vector 
%  Initialize  Yaw  rate  vector 

%  Initialize  X  position  vector ( including  side  slip) 

%  Initialize  forward  velocity  vector (including  side  sli 

%  Initialize  Y  position  vector ( including  side  slip) 

%  Initialize  side  slip  velocity  vector ( including  side  s 


%  Define  storage  vectors  for  Side  Slip  Model 
% 

zeros { 1 , kmaxl ) ;  %  Initialize  side  slip  veloctiy  vector 

V  dot  -  zeros  (1 , kmaxl) ;  %  Initialize  side  slip  acceleration  vector 

X  ^  Initialize  X  position  vector  ( including  side  slip) 

X2dot  =  zeros  (1, kmaxl);  %  Initialize  forward  velocity  vector ( including  side  sli 

^  ^  Initialize  Y  position  vector  { including  side  slip) 

Y2dot  =  zeros  (1,  kmaxl);  %  initialize  side  slip  velocity  vector  ( including  side  s 

% 

%  Define  storage  vectors  for  Nomoto's  Model  w/correction  factor 

zeros ( 1 , kmaxl ) ;  %  Initialize  angular  velocity  vector 

r3dot  -  zeros  (1, kmaxl) ;  %  Initialize  angular  acceleration  vector 

psi3  -  zeros  (1 , kmaxl ) ;  %  initialize  heading  angle  vector 

psi3dot  -  zeros ( 1 , kmaxl ) ;  %  Initialize  Yaw  rate  vector 

^  Initialize  X  position  vector  ( including  side  slip) 
X3dot  -  zeros  (1 ,  kmaxl)  ;  %  Initialize  forward  velocity  vector  (including  side  sli 

Y3  =  zeros { 1 , kmaxl )  ; 

Y3dot  =  zeros { 1 , kmaxl ) ; 
lip) 

Timel  =  zeros { 1 , kmaxl); 


%  Initialize  Y  position  vector ( including  side  slip) 

%  Initialize  side  slip  velocity  vector ( including  side  s 

%  Initialize  time  vector 


for  i  =  1 : kmaxl- 1; 

% 

%  Calculation  loop  for  Nomoto's  Model 
% 

rldot(i)  =  (K*d  -  rl{i))/T; 
rl(i+l)=  rl{i)  +  rldot(i)*dt; 
psildot (i)  =  rl (i) ; 
psil(i+l)=  psil{i)  +  psildot (i) *dt; 

XLdot (i)  =  cos (psil { i ))  ; 

XL(i+l)  =  XL(i)  +  XLdot (i)^dt; 

YLdot (i)  =  sin (psil { i  ))  ; 

YL(i+l)  =  YL(i)  +  YLdot (i)*dt; 

% 

%  Calculation  loop  for  Side  Slip  Model 
% 

v2dot(i)  =  (Kv*d  -  v2(i))/Tv; 
v2(i+l)  =  v2(i)  +  v2dot ( i ) *dt ; 

X2dot(i)  =  cos  (ps'il  (i)  )  -  v2  { i )  *sin  (psil  ( i )  )  ; 

X2{i  +  1)  =  X2(i)  +  X2dot  (i)  “^dt; 

Y2dot{i)  =  sin(psil (i) ) +  v2 ( i ) *cos (psil ( i ) ) ; 

Y2{i+1)  =  Y2(i)  +  Y2dot(i)*dt; 

% 

%  Calculation  loop  for  Nomoto's  Model  w/correction  factor 
% 

R=0.5424*d+1.2560; 


r3dot(i)  =  (R*K*d  -  r3(i))/T; 
r3(i+l)=  r3(i)  +  r3dot(i)*dt; 
psi3dot (i)  =  r3 (i) ; 
psi3.(i  +  l)=  psi3(i)  +  ps  i3dot  (  i  )  *dt  ; 
X3dot{i)  =  cos (psi3 ( i ) ) ; 

X3{i+1)  =  X3(i)  +  X3dot(i)*dt; 
Y3dot(i)  =  sin (psi3 ( i )  )  ; 

Y3(i+1)  =  Y3(i)  +  Y3dot{i)*dt; 

Timel (i+1)  =  Timel (i)  +  dt; 
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% 

end 

% 

%  NONLINEAR  MODEL 
% 

%  Define  hydrodynamic  coefficients 
% 

Yb  =  -0.5*PI*k  +  1.4*s*Cb; 

YR  =  ml  -  1.5*s*Cb; 

Ybb  =  1 .25* (k/s) * (1-Cb)  +  0.5; 

Yrr  =  0 . 17* (k/s) *Cb  -  0.07; 

Ybrr  =  2.97* (k/s) * (1-Cb) ; 

Ybbr  =  0 .75* (k/s) *Cb  -  0.65; 

Nb  =  -k; 

NR  =  -0.54*k  +  k^2; 

Nbb  =  -0.48* (k/s) * (1-Cb)  +  .066; 

Nrr  =  0.5*s*Cb  -  .09; 

Nbrr  =  -0 . 25* (k/s ) *Cb  +  0.05; 

Nbbr  =  -57.5*  (s*Cb)"'2  +  18.4*s*Cb  -  1.6; 

Xudot  =  m-ml ; 

Nrdot  =  Izz-I; 

Yvdot  =  m-m2 ; 

Yrdot  =  0.0; 

Nvdot  =  0.0; 

Xbr  =  0.0; 

Xuu  =  0.0; 

U  =  1.0; 

UR  =  U;  %  effective  rudder  inflow  speed 


%  Define  storage  vectors  in  advance  to  optimize  performance 
% 

kmax  =  (time/dt)+l;  %  Number 

V  =  zeros ( 1 , kmax) ;  %  Initial 

r  =  zeros ( 1 , kmax) ;  %  Initial 

rdot  =  zeros  ( 1 ,  )<max)  ;  %  Initial 

b  =  zeros ( 1 , kmax)  ;  %  Initial 

bdot  =  zeros ( 1 , kmax)  ;  %  Initial 

X  =  zeros ( 1 , kmax)  ;  %  Initial 

Xdot  =  zeros { 1 , kmax)  ;  %  Initial 

Y  =  zeros { 1 , kmax)  ;  %  Initial 

Ydot  =  zeros ( 1 , kmax)  ;  %  Initial 

psi  =  zeros { 1 , kmax)  ;  %  Initial 

psidot  =  zeros ( 1 , kmax) ;  %  Initial 

Time  =  zeros  { 1 ,  )anax)  ;  %  Initial 

% 

%  Define  nonlinear  equations  [calculation  loop] 

% 

for  i  =  l;kmax-l; 


of  data  points 

ize  side  slip  veloctiy  vector 

ize  angular  velocity  vector 

ize  angular  acceleration  vector 

ize  drift  angle  vector 

ize  drift  acceleration  vector 

ize  X  position  vector 

ize  forward  velocity  vector 

ize  Y  position  vector 

ize  side  slip  velocity  vector 

ize  heading  angle  vector 

ize  yaw  rate  vector 

ize  time  vector 


- (AR/ (L*D) ) *Cn*UR^2* 
- { 1-tr ) *Fn*sin(d) ; 

- { 1+aH) *Fn*cos (d) ; 
-(xR  +  aH*xH) *Fn*cos 


Fn 
Xr 
Yr 
Nr 
% 

Xh  =  Xbr *r ( i ) *sin (b { i )  ) 
Yh  =  b(i)*Yb  +  r ( i ) *YR  + 
) * (r (i) ) ^2*Ybrr; 

Nh  =  b(i)*Nb  +  r(i)*NR  + 
) * (r (i) ) "2*Nbrr; 

% 


sin(d  -  gamma* (b ( i ) - ( 2 *xR*r ( i ) ) ) ) ;  %  Force  acting  on  ruddt 


:d)  ; 

-  Xuu  *  (cos (b ( i)  )  ) ^2  ; 

r ( i) * (abs (r (i) ) ) *Yrr  +  b ( i ) * (abs (b ( i ) ) ) *Ybb  +  (b ( i ) ) ^2*r ( i ) *Ybbr 
r  (i)  *  (abs  (r  (i)  )  )  *Nrr  +  b ( i )  *  (abs  (b  (  i )  )  )  *Nbb  +  (b  (  i  )  ) '"2*r  (  i  )  *Nbbr 


Yp  =  0.0; 

Np  =  0.0; 

Jp  =  (U*cos (b ( i ) ) * (1-wp) ) / (n*Dp) ; 

cl=0.52; 

c2=-0.4861; 

c3  =  0. 01212; 

Kt  =  cl  +  c2*Jp  +  c3*Jp^2; 

Xp  =  (Ctp*  (1-tp)  *n"2*Dp"'4*Kt)  /  (0.5*L*D*U^2)  ; 

% 

XI  =  Xh  +  Xr  +  Xp; 

Y1  =  Yh  +  Yr  +  Yp; 

N1  =  Nh  +  Nr  +  Np; 

% 

fl  =  (m- Yvdot )*( Izz-Nrdot)  -  {Nvdot-m*Xg ) * ( Yrdot-m*Xg ) ; 
f 2  =  Y1 -  m*r ( i ) ; 
f3  =  Nl-  m*Xg*r(i); 

% 

rdot(i)  =  [  (m- Yvdot) *f3  +  (Nvdot-m*Xg) * f 2 ] / f 1  ; 


+  b( 
+  b( 
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r(i+l)  =  r(i)  +  rdot(i)*dt; 

bdot(i)  =  [ (Izz-Nrdot) *f2  +  (Yrdot-m*Xg) * f3 ] /f 1 ; 
b{i+l)=  b(i)  +  bdot{i)*dt; 

V  { i )  =  U* tan (b ( i ) ) ; 

psidot ( i)  =  r  (i) ; 

psi(i+l)  =  psi(i)  +  psidot (i) *dt; 

Xdot{i)  =  cos{psi(i))  -  v(i) *sin(psi (i)  )  ; 

X(i+1)  =X(i)  +  Xdot{i)*dt; 

Ydot(i)  =  sin(psi(i))  +  v  ( i ) *cos (psi ( i )  )  ; 

Y(i+1)  =  Y(i)  +  Ydot{i)*dt; 

Time(i+1)  =  Time(i)  +  dt; 

end 

% 

%  Plot  results 
% 

figure ( 1 ) 

plot {Y,X,  ' , YL,XL.  ' - ' , Y3 ,X3,  '  .  '  ) , grid; 
xlabel ( ^Transfer' ) , ylabel ( 'Advance' )  ; 
gtext('-“  Nonlinear  Model'); 
gtext ( '  .  Nomoto  model  w/correction '  )  ; 
gtext ( ' -  Nomoto  Model ' ) ; 

% 

figure (2) 

plot (Time (1 : 101),r(l:101),'--' ,Timel (1 : 101 ) , rl ( 1 : 101)  ,  '  - ' ,Timel (1 : 101) ,r3  (1:101) , '  . ' )  .grid 

xlabel ( ' Time (sec ) ' ) , ylabel ( ' Yaw  Rate ' ) 

gtext (' —  Nonlinear  Model'); 

gtext ( ' .  Nomoto  model  w/correction'  )  ; 

gtext (' -  Nomoto  Model'); 

% 
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